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We present preliminary results of our analysis of the Fermi -LAT data from the direction of NGC 4151. We 
find a new 7 -ray source with a statistical significance cr > 5, shifted by 0.5° from the position of NGC 4151. 
Apparently, the source was bright only during a 1.5-year period between December 2011 and June 2013 and it 
strongly contaminated the signal from NGC 4151. Therefore, we neglect this period in our analysis. We find 
two additional, persistent 7 -ray sources with high cr, shifted from NGC 4151 by ~ 1.5° and 5°, whose presence 
has been recently confirmed in the Third Fermi Catalog. After subtracting the above sources, we still see a weak 
residual, with cr < 3, at the position of NGC 4151. We derive an upper limit (UL) for the 7 -ray flux from NGC 
4151 and we compare it with predictions of the ADAF model which can explain the X-ray observations of this 
object. We find that the Fermi UL strongly constrains non-thermal acceleration processes in hot flows as well 
as the values of some crucial parameters. Here we present the comparison with the hot flow models in which 
heating of electrons is dominated by Coulomb interactions with hot protons. In such a version of the model, 
the 7 -ray UL, combined with the X-ray data, constrains the energy content in the non-thermal component of 
proton distribution to at most a few per cent, rules out a weak (sub-equipartition) magnetic field and favors a 
rapid rotation of the supermassive black hole. 


1. INTRODUCTION 

Low-luminosity AGNs, with the luminosities below 
~ O.OlLEdd, are likely to be powered by optically 
thin, hot accretion flows (a.k.a. ADAFs, see e.g. [12]). 
The two-temperature structure is a key property of 
ADAFs, as such flows are supported by the proton 
pressure. In their innermost parts, the hot protons 
have energies above the threshold for pion production. 
As estimated e.g. in [6,8], the decay of pions leads to 
substantial fluxes of 7-rays, which may be probed in 
nearby AGNs at the current sensitivity of Fermi -LAT 
surveys. The Fermi -LAT data for radio-quiet AGNs 
were analyzed in [10] (two years of data) and [2] (three 
years) and the derived upper limits are already quite 
stringent compared to expectation. In [ll]we revisit 
the issue of searching the signatures of hadronic emis¬ 
sion from hot flows and the related implications for 
hot-flow models. We perform the detailed analysis of 
nearby, low-luminosity Seyfert galaxies using over 6 
years of Fermi -LAT data and we compare the results 
with the model predictions for a complete range of the 
model parameters. In this contribution we report our 
preliminary results for one of the best-studied AGNs, 
NGC 4151. 

NGC 4151 is one the X-ray brightest AGNs, with 
the bolometric L ~ O-Ol+Edd, showing no signa¬ 
tures of a relativistically distorted reflection compo¬ 
nent (constraints on the width of Fe Ka line imply 
the lack of an optically thick material within at least 
the innermost ~ 100I?g) as well as showing a spectral 
similarity to black hole binaries in their hard states 
(also most likely powered by hot flows), see e.g. [4]. 
All these properties make it a relevant objects for test¬ 


ing the hot flow scenario. Below we use the black hole 
mass M = 3.8 x 1O 7 M 0 from the stellar dynamical 
mass measurement [9]. 

2. LAT DATA ANALYSIS 

We analyzed the data from the direction of NGC 
4151, comprising 6.4 years of Fermi-LAT observations 
carried out between 2008 August 4 and 2015 January 
10. Events we selected from a region with the ra¬ 
dius of 10° centered on the position of NGC 4151. 
We performed the unbinned likelihood analysis us¬ 
ing the v9r33p0 Fermi Science Tools with CALDB 
instrument response functions. We used the standard 
templates for the Galactic (gll_iem_v05_revl. f its) 
and the isotropic (iso_source_v05jrevl.txt) back¬ 
grounds. In our initial model of the region we took 
into account only the sources from the Second Fermi 
Catalog (2FGL) [7], i.e. our model included the same 
sources as those used in [2]. 2FGL ,11209.6+4121 
(marked by the red circle in Fig. la) is the 2FGL 
source closest to NGC 4151, with the distance of ~ 2°. 

Fig. la shows the TS map of the region, built af¬ 
ter subtracting the 2FGL sources. The map reveals 
residual structures indicating the presence of addi¬ 
tional point-like sources, marked by the green cir¬ 
cles. For each of these objects we use the gtlike and 
gtfindsrc tools to find its significance, best-fit posi¬ 
tion and spectral parameters; the results are given in 
TableQ] Sources SI (beyond the map in Fig. 1) and S3 
have been recently reported in 3FGL [1], with param¬ 
eters very similar to these estimated in our analysis. 
S2 is not reported in 3FGL, however, this source is 
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Figure 1: Test-statistic maps of the 5° x 5° region around the position of NGC 4151 (yellow circle) with the pixel size of 
0.125° (left) and 0.15° (right). The left map was generated from the data from the total 6.4 years and only the sources 
from 2FGL were subtracted from this map (the only 2FGL source in this region of the sky is marked by the red circle); 
the green circles indicate new sources revealed in our analysis. The right map was generated from the data taken before 
December 2011 and after June 2013 and the sources from 3FGL (red circles) were subtracted from this map. 


Table I New sources introduced in the model of the 
region around NGC 4151, not reported in 2FGL. (1) 
Source (see Fig. 1), (2) 3FGL name, (3) gtlike TS 
values, (4) and (5) gtf indsrc coordinates. Results for S2 
were obtained using the data taken between December 
2011 and June 2013; for the remaining sources the total 
data set for 6.4 years was used. 


(i) 

(2) 3FGL name 

(3) TS 

(4) RA 

(5) DEC 

si 

3FGL J1220.2+3434 

191 

185.06 

34.57 

S2 

not reported 

30 

182.86 

38.95 

S3 

3FGL J1203.2+3847 

32 

180.81 

38.79 

S4 

not reported 

22 

184.90 

36.93 


critical for the analysis of a signal from NGC 4151, 
as its distance of ~ 0.5° is comparable (or smaller 
below ~ 1 GeV) to the LAT point spread function. 
Therefore, we check properties of this source in more 
details. 

By using gtsrcprob we find that the position of 
S2 is determined mostly by 4 photons with energies 
between 10 and 20 GeV which arrived from the same 
direction (within 10 arcrnin) between December 2011 
and June 2013. At lower energies, neglecting the four 
events with E > 10 GeV, we also see the signature of 
increased activity of S2 during that period, in the form 
of an extended residual covering the nominal positions 
of NGC 4151 and S2. 

In the TS map built for the data neglecting the 
above period (see Fig. lb) we do not see a strong 


signal at the S2 position, we therefore conclude that 
S2 strongly dominated the emission from the region 
around NGC 4151 only during the 1.5 year out of the 
total 6.4 considered years. Note that S2 is not re¬ 
ported in 3FGL which includes sources detected with 
TS > 25 using the data taken during the four years 
up to 2012 July (i.e. covering only ~ 30% of the time 
of the increased activity of S2). 

For our further analysis we neglect the data taken 
during the 1.5 year when S2 was bright. We subtract 
the sources reported in 3FGL and we get the TS map 
shown Fig. lb. The map shows a weak residual, which 
can be fully compensated for by adding the source at 
the nominal position of NGC 4151 and gtlike gives 
TS ~ 8 for such a source. At its very low statistical 
significance, it is not possible to assess whether it rep¬ 
resents a background fluctuation or an actual emission 
from the studied object; it may also contain some con¬ 
tribution from emission of S2 in its lower luminosity 
states. 

We then derive the 95% confidence level upper limit 
(UL) for the integrated photon flux from NGC 4151 
neglecting the data between December 2011 and June 
2013. The pion decay spectra can be approximated by 
a simple power-law only in limited energy ranges (see 
Fig. 2 below), therefore, we assume relevant values of 
the photon index, T, and find the UL in the 0.3-1 
GeV range to compare with the 7r°-decay spectra for 
the thermal distribution of protons, and in the 1-10 
GeV to compare with the model assuming a power- 
law distribution of protons. The results are given in 
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Figure 2: Spectra received by a distant observer from a hot flow around the black hole with M = 4 x 10 7 Mq for 
a = 0.95, p = 9, 5 = 10 -3 and m = 0.3. The black solid line shows the component formed by the cooling of electrons 
(mostly thermal Comptonization). The blue and red solid lines show the component formed by 7 r°-decay for thermal 
and power-law (with s = 2 . 1 ) distribution of protons, respectively; the thin dashed lines show the corresponding 
spectra of the hadronic component in the rest frame of the flow (the spectra shown by solid lines are affected by the 
GR transfer and 77 absorption). 


Table II 95% C.L. upper limits for the photon flux (3) 
and 7 -ray luminosity (4) for the energy ranges given in 
column ( 1 ) and the assumed value of 7 given in 
column ( 2 ). 


(1) Energy 
range [GeV] 

(2) Photon 
index 

(3) UL: F 
[phot /cm 2 /s] 

(4) UL: L 7 
[erg/s] 

0.3 - 1 

4 

6 x 10 _1 ° 

6.5 x 10 39 

1 - 10 

2.1 

1.5 x 10 _1 ° 

3.1 x 10 39 

1 - 10 

2.7 

1.3 x 10 ~ 10 

1.9 x 10 39 


Table HI 


3. GAMMA-RAY EMISSION FROM HOT 
FLOWS AND COMPARISON WITH NGC 
4151 

In [5,6] we present a precise model for emission 
from hot flows, taking into account the relevant lep- 
tonic and hadronic processes, and using ( 1 ) a fully 
general relativistic (GR) description of both the ra¬ 
diative and hydrodynamic processes; ( 2 ) an exact, 
Monte Carlo computation of global Comptonization; 
(3) seed photons input from nonthermal synchrotron 
emission of 7 r -decay electrons; (4) an exact computa¬ 
tion of the absorption of 7 -ray photons in the radiation 
field of the flow. The model is parametrized by: the 
black hole mass, M, the dimensionless accretion rate, 
m = Mc 2 /LEdd, the ratio of gas to magnetic pressure, 
/3, the fraction of the dissipated energy which heats 
directly electrons through MHD processes, <5, and the 
spin parameter, a = J/(cR s M), where J is the black 
hole angular momentum and R s = GM/c 2 . We take 
into account models with thermal and with power-law 
distrubutions of protons; in the latter the power-law 
index, s, is also a free parameter. 


Here we focus on models with small 5, i.e. with 
electrons heated by Coulomb interactions. We briefly 
summarize properties crucial for our final conclusions; 
Fig.[2]shows example spectra of radiation produced in 
a hot flow by thermal Comptonization and by 7 r° de¬ 
cay. 

(i) The nonthermal synchrotron radiation from 
decay electrons gives the dominating input of seed 
photons for Comptonization and it allows to reconcile 
the hot-flow model with the AGN X-ray data. It also 
provides an attractive explanation of spectral differ¬ 
ences between AGNs and black-hole transients within 
the same physical model, see [5]. 

(ii) For bolometric L ~ (0.001 — 0.01)LEdd, the size of 
the 7 -ray photosphere (inside which the flow is opaque 
to 7 -rays) equals several R g . As a result, for models 
assuming thermal protons, the 7 -ray flux detected by 
a distant observer is reduced by several orders of mag¬ 
nitude, because the 7 -rays are produced mostly inside 
the photosphere. We note that [ 8 ] assessed compa¬ 
rable X-ray and 7 -ray fluxes from flows surrounding 
rapidly rotating black holes. However, they neglected 
the GR transfer and 77 absorption; taking into ac¬ 
count these effects we get, for thermal protons, the 
7 -ray fluxes smaller by ~ 3 (large /?) to ~ 5 (small /3) 
orders of magnitude than the X-ray flux. 

(iii) In models assuming a thermal distribution of pro¬ 
tons, the 7 -ray flux is extremely sensitive to the value 
of /3. In flows with smaller /3 (larger B), a larger frac¬ 
tion of the accretion power is used to build up the 
magnetic field strength; therefore, the energy heating 
the particles, and hence the proton temperature, is 
smaller. As for thermal protons the 7 -ray luminosity, 
L 7 , is extremely sensitive to the proton temperature, 
the above effect leads to the difference by 2 - 3 orders 
of magnitude between L 1 for a strong (in equipartition 
with gas) and weak magnetic field. 

(iv) The proton temperature increases with a and 
hence the 7 -ray emissivities strongly depends on a for 
thermal protons. However, the largest difference be- 
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Figure 3: The 7 -ray (1-10 GeV) luminosity Eddington 
ratio as a function of the X-ray (2-10 ken) luminosity 
Eddington ratio. The squares and circles show the hot 
flow model predictions for the nonthermal proton 
distribution with s = 2.1. The Fermi UL was obtained for 
the assumed Y = 2.1 in the 1-10 GeV range (see Table 2) 
and the average X-ray luminosity was estimated from the 
Swift-BAT data. The red circles are for a = 0.998, /3 = 1, 
m = 0.1 and 0.3; the blue squares are for a = 0.95, /3 = 9, 
m = 0.3 and 0.8. 



Figure 4: The same as in Fig. 3 but the model points are 
for the thermal distribution of protons and the 7 -ray 
Eddington ratio is determined for the 0.3-1 GeV range; 

T = 4 was assumed for the Fermi UL. 


tween the emissivities occurs within the photosphere 
and, therefore, the dependence of the observed L T on 
a is reduced by 77 absorption. On the other hand, 
a sufficiently strong input of seed photons from the 
emission of 7 r -decay electrons requires either a rapid 
rotation of the black hole or a significant content of 
nonthermal protons. 

Using the publicly available Swift-BAT data [3] we 
find the average X-ray luminosity of NGC 4151 be¬ 
tween 2009 and 2014 (we find that it corresponds to 
the intermediate state as defined in [4]), which then 
allows us to compare the L T predicted by the model 
with the Fermi UL. The results are shown in Figs 3 
and 4. For the power-law distribution of protons (Fig. 
3) the UL is over an order of magnitude lower than the 
model prediction, which constrains the energy content 
in the nonthermal component of the proton distribu¬ 


tion to at most a few per cent. 

For the thermal distribution of protons (Fig. 4), the 
Fermi UL is sufficiently low to exlude models with 
large /3, which, taking into account the above, rules 
out any version of the model with a weak magnetic 
field. For an equipartition value of /3 ~ 1 the predicted 
flux is below the UL value. 

Then, the Fermi data favor a strongly magnetized 
plasma with a weak content of nonthermal protons. 
For such a case, a high spin value is required for a suf¬ 
ficiently strong flux of seed photons from nonthermal 
emission of pion-decay electrons. 

4. SUMMARY 

We thoroughly analyzed the 7 -ray data from a re¬ 
gion around NGC 4151, which led us to identification 
of new 7 -ray sources. After subtracting their contri¬ 
bution, we get a weak residual signal at the position of 
NGC 4151. At its low statistical significance, a < 3, 
it is not possible to assess its nature. 

Comparison of the derived upper limits with the 
model predictions allows to constrain several crucial 
quantities which illustrates the potential of Fermi 
measurements in probing the properties of flows pow¬ 
ering AGNs at low luminosities. 
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